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ABSTRACT 


Presented  is  a  method  to  estimate  the  element  gains  and  phases 
of  a  sensor  array  when  interfering  signals,  including  specular 
multipath,  are  present.  A  movable  calibration  source  is  used.  The 
element  gains  are  derived  from  the  relative  output  powers  at  the 
array  elements.  The  element  phases  are  obtained  by  transforming  the 
array  snapshots  to  make  the  source  direction  appear  fixed  and  then 
maximizing  a  function  constructed  with  the  transformed  snapshots. 

Computer  simulations  produce  highly  accurate  estimates  of 
element  gains  and  phases  when  the  signal  directions  are  known  and 
the  signal -to-noise  ratio  is  reasonably  high.  When  the  directions 
are  not  known,  the  set  of  element  phase  estimates  contains  an 
additional  component  which  produces  a  bias  in  direction  estimates. 


RfiSUM^ 

Nous  prdsentons  une  methods  pour  estimer  les  gains  et  les 
phases  des  414ments  d'une  antenne  rdseau  en  presence  de  brouillage 
incluant  la  reflexion  sp4culaire.  Une  source  d'Stalonnage  amovible 
est  utilisee.  Les  gains  des  414ments  sont  obtenus  par  les 
puissances  relatives  des  sorties  des  §16ments  du  reseau.  Les  phases 
sont  obtenues  en  maximisant  une  fonction  construite  avec  les 
vecteurs  de  donn4es  transformes  de  telle  sorte  que  la  direction  de 
la  source  apparaisse  fixe. 

Quand  la  direction  des  signaux  est  connue  et  que  le  rapport 
Signal -sur- Bruit  est  raisonnablement  Sieve,  nous  obtenons  des 
estimSs  trSs  prScis  des  gains  et  des  phases  des  SlSments  par  des 
simulations  sur  ordinateur.  Par  contre,  quand  les  directions  ne 
sont  pas  connues,  1' ensemble  des  estimSs  des  phases  contient  une 
composante  supplSmentaire  qui  produit  un  biais  dans  I'estimS  des 
directions. 


EXECUTIVE  SUMMARY 


Accurate  calibration  of  a  sensor  array  is  essential  when  this 
array  is  used  in  direction- finding  with  modern  high-resolution 
signal  processing  methods  such  as  the  MUSIC  method. 

This  report  discusses  a  new  calibration  technique  that 
estimates  the  element  gains  and  phases  of  an  array  when  interfering 
signals,  including  specular  multipath,  are  present.  The  technique 
uses  a  movable  pilot  source  to  generate  the  calibration  data. 
Initially,  the  element  gains  are  estimated  from  the  output  powers 
at  the  array  elements.  Next,  a  transformation  is  used  to  equalize 
the  element  gains  and  simultaneously  fix  the  source  direction  in 
the  data.  Finally,  the  element  phases  are  estimated  by  maximizing 
an  objective  function  constructed  with  the  transformed  data. 

Computer  simulations  show  that  the  estimates  of  the  element 
gains  and  phases  are  highly  accurate  if  the  calibration  source 
directions  are  known  and  the  signal-to-noise  ratio  at  the  array 
elements  is  30  dB  or  higher. 
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1 . 0  INTRODUCTION 


Accurate  calibration  of  a  sensor  array  is  essential  if  the 
outputs  of  this  array  are  to  be  used  in  direction- finding  with 
high-resolution  methods  such  as  the  MUSIC  method  [1] .  Otherwise, 
the  ed^ility  of  these  methods  to  resolve  signal  directions  would  be 
degraded.  A  discussion  of  how  sensor  errors  degrade  the  performance 
of  the  MUSIC  method,  for  example,  is  given  by  Friedlander  in  [2] . 

Many  techniques  to  calibrate  a  narrowband  sensor  array  in  the 
absence  of  specular  multipath  interference  have  been  developed. 
They  include  those  by  Paulraj  and  Kailath  [3]  ,  Friedlander  and 
Weiss  [4,5],  Weiss,  Willsky,  and  Levy  [6],  Brown,  McClellan,  and 
Holder  [7,8],  and  Pierre  and  Kaveh  [9].  In  [3]  to  [6]  the  authors 
are  concerned  with  the  estimation  of  element  gains  and  phases  when 
mutual  coupling  among  the  array  elements  can  be  ignored.  References 
[7]  to  [9]  include  the  effects  of  mutual  coupling. 

The  condition  that  specular  multipath  interference  be  absent 
is  frequently  not  satisfied  when  the  array  has  to  be  calibrated  at 
the  site  of  operation.  This  interference  is  unavoidable,  for 
example,  in  the  special  case  of  an  array  antenna  used  to  track  low- 
flying  targets  over  water,  a  high-frequency  array  used  in 
direction- finding,  and  a  sonar  array  for  submarine  detection. 

This  paper  is  concerned  with  the  estimation  of  element  gains 
and  phases  when  specular  multipath  and  other  types  of  interference 
are  present.  A  movable  calibration  source  is  used.  Motion  of  the 
sensor  array  is  allowed.  Computer  simulations  show  that  the  element 
gain  and  phase  estimates  are  highly  accurate  if  the  calibration 
source  directions  are  known  and  the  signal -to-noise  ratio  (SNR)  at 
the  array  elements  is  reasonably  high. 

The  theory  for  this  array  calibration  method  are  derived  in 
Section  2 .  The  procedure  is  described  in  Section  3 .  Some  results  of 
performcuice  evaluations  are  presented  in  Section  4 .  The  conclusions 
are  given  in  Section  5. 
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2 . 0  THEORY 


This  estimation  method  is  based  on  the  use  of  a  reasonably 
strong  movable  calibration  source.  The  source  motion  is  exploited 
to 

(a)  Separate  the  estimation  of  element  gains  from  the 
estimation  of  element  phases, 

(b)  Decorrelate  the  direct  and  indirect  calibration  signals 
at  the  array  elements,  and 

(c)  Make  each  interfering  signal  source  behave  like  a  source 
spread  out  over  a  range  of  directions. 

Some  examples  of  source  motion  will  be  given  later. 

In  the  notation  used,  N  denotes  the  number  of  array  elements; 
a,  the  vector  (a^, aj,  •  •  • , a^) diag{b},  the  diagonal  matrix 
diag{bi,b2,  •  •  •  ,bH} ;  {gi,g2/ •  •  *  .gw} »  the  element  gains;  {0i, 02#  •  *  •  / 0n}  / 
the  element  phases;  [*  1  ,  the  absolute  value;  and  E{*  },  the  expected 
value.  Superscripts  T  and  H  represent  the  transpose  and  conjugate- 
transpose,  respectively.  The  uniqueness  conditions  for  the  element 
gains  and  phases  are  defined  as 

g2  +  g2+..-+g2  =  N  (1) 


and 


(f)i  +  <p2  =  0  .  (2) 

The  remainder  of  this  section  is  organized  as  follows . 
Initially,  the  signal  model  and  the  estimation  of  element  gains  are 
described.  Next,  it  is  shown  that  there  is  a  transformation  which 
equalizes  the  element  gains  in  the  calibration  data  and 
simultaneously  make  the  calibration  source  direction  appear  fixed. 
Finally,  the  estimation  of  element  phases  is  discussed. 

2 . 1  Signal  Model 

The  array  elements  are  isotropic.  Mutual  coupling  is  absent. 
Given  that  the  calibration  source  direction  is  9^,  a  snapshot  of 
the  array  element  output  is 

S  =  a  +  interference  +  noise  ,  (3) 

where  a  is  the  signal  amplitude  and  b(g,^,  is  the  array  manifold 
vector  for  6^.  An  array  manifold  vector  for  9  has  the  decomposition 

b(g,<fe,fl)  =  G(g)  T(^)  a(9)  ,  (4) 
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where 


fi(3)  =  diag{g}  , 
l(^)  =  diag{x{^)}  , 


(5) 


(6) 


vector  i(^)  is  defined  as 

Tn(±)  -  exp[j0J  ,  (7) 

and  ^{6)  is  the  array  steering  vector.  For  the  sake  of  convenience, 
it  will  be  assumed  that  the  components  of  a(0)  have  unit 
magnitudes . 

2 . 2  Estimation  of  Element  Gains 

By  design,  the  direct  signal  from  the  calibration  source  is 
very  strong  compared  with  the  signals  from  other  sources.  Specular 
multipath,  if  present,  is  effectively  uncorrelated  with  the  direct 
signal  at  the  array  elements.  The  element  gains  g^  to  g^  can 
therefore  be  derived  from  the  average  output  powers  at  the  array 
element.  This  calculation  is  straightforward  and  is  given  in  the 
estimation  procedure  in  Section  III. 

The  direct  and  specular  multipath  signals  are  not  correlated 
because  their  phase  difference  at  each  array  element  changes  when 
the  calibration  source  direction  changes.  Careful  planning  can  make 
this  difference  change  by  thousands  of  degrees  in  the  data. 

2.3  Transformation  for  Calibration  Data 

The  array  manifold  vector  for  the  calibration  source  direction 
is  given  by  (4)  with  8=6^.  A  simple  rearrangement  yields 

=  S(a)  A(flc)  i(^)  /  (8) 


where 


a(0c)  =  diag{a(flc)}-  (9) 

This  rearrangement  makes  the  phases  of  a(0c)  behave  like  element 
phases.  It  also  makes  t{^)  behave  like  an  array  steering  vector. 

Let  a  transformation  matrix,  dependent  on  8^,  be  defined  as 

=  A'M^c)  S'Mg)  .  (10) 

Applied  to  the  array  manifold  vector  given  by  (4) ,  it  gives 

r(flc)  b(g,<fe,(9)  =  A'Mflc)  A(0)  T(^) 

=  zCife)  /  if  ^  =  •  (11) 
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Because  ^(^c)  and  ^{d)  are  diagonal  and  the  diagonal  elements  have 
unit  magnitudes,  the  transformed  manifold  vector  can  be  treated  as 
a  vector  with  unit  element  gains,  zero  element  phases,  and  a  new 
steering  direction  determined  by  the  right  hand  side  of  (11) .  For 
the  new  direction  is  determined  by  the  phases  of  and  is 

independent  of  8^. 

Transformation  matrix  £(dc)  is  designed  for  the  calibration 
data  in  which  each  array  snapshot  has  a  different  8^.  It  spreads 
out  the  power  of  each  interfei'ing  source  over  a  range  of  directions 
and  simultaneously  make  the  calibration  signal  behave  like  the 
strongest  localized  signal  in  the  data. 

2.4  Estimation  of  Element  Phases 

Matrix  E(dc)  converts  an  array  snapshot,  given  by  (3) ,  to 

y  =  £(^c)  X 

=  a  X(^)  +  z  ,  (12) 

where  3,  is  the  contribution  from  interference  and  noise.  Because 
the  calibration  signal  is  uncorrelated  with  the  interfering 
signals,  the  correlation  matrix  constructed  with  y  is 

Sy  =  E{y 

=  lap  X(^)  i“(i)  +  S*  (13) 

where  Rj=E{zz"} .  In  Ry,  the  element  gains  are  unity,  the  element 
phases  are  zero,  and  the  array  steering  vector  for  calibration 
source  direction  is  x(i) • 

Let  the  Capon  method  [10,11]  be  used  to  estimate  the  array 
output  power  for  the  source  direction  in  Ry.  If  ^  is  an  estimate  of 
^  in  this  calculation,  the  power  is 

1 

<P(^,Ey)  =  -  .  (14) 

F(^,Ry,i) 

where 

F(^,Ry,^)  =  x“(<6a)  Sy  X(^)  •  (15) 

Because  Ry  is  given  by  (13) ,  application  of  the  matrix  inversion 
lemma  leads  to 

lap  Rz^x(^)  x"(^)  Sz^ 

R-^i  =  r;1 - .  (16) 

1  +  lofp 
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and 


<P(i,Sy.) 


1  +  Ia|"  F(^,R„(ji) 
F(i,R„i)  +  lap  G{^,S*,i) 


(17) 


with 

G(^,E.,i)  -  F(^.fi.,^)  F(i,E*,i)  -  |f(<^..S*,s^)  M 

=  IsM  1^1"  -  la”  aM  .  (18) 

The  second  line  is  obtained  by  noting  that  R;^  is  the  inverse  of  a 
non- singular  correlation  matrix  when  noise  is  present  and  that  it 
has  the  decomposition 

V”  .  (19) 

where  V  is  a  non- singular  matrix  also.  Vectors  g  and  ^  are  defined 
as 

a  =  V  iv^)  (20) 

^  =  V  T(i)  .  (21) 

It  follows  from  (2),  (7),  and  Schwarz's  inequality  that 

G(<fe,R,,i)  2  0  ,  (22) 

The  equality  sign  holds  if  and  only  if 

The  above  results  indicate  that  an  estimate  of  ^  can  be 
obtained  by  maximizing  (P(i,Ry)  with  respect  to  provided  that  the 
rate  at  which  changes  in  the  neighbourhood  of  ^  is  small 

compared  with  that  of  lo(pG(^,Rj,^)  .  Increasing  |al  yields  a  more 
accurate  estimate  of  Making  F(^,Rj,^)  independent  of  ^  yields 
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3.0  ESTIMATION  PROCEDURE 


There  are  two  variations  of  the  calibration  procedure: 
calibration  source  directions  known,  and  calibration  source 
directions  not  known.  Both  involve  the  maximization  of  an  objective 
function  with  respect  to  the  estimates  of  element  phases .  Some 
general  properties  of  this  fvinction  are  also  discussed. 

3 . 1  Algorithm  1 :  Source  Directions  Known 

Given  a  set  of  M  array  snapshots  {iSi/Xj,  •  •  • , and 
corresponding  calibration  source  directions  { ^i/ ^2/ ' ' " / ®m}  /  the 
estimation  of  array  element  gains  and  phases  is  as  follows: 

Step  1. 

Estimate  the  average  output  powers  at  the  array  elements  as 
IxinM  +  IxjnM  +  •  •  •  +  Ix^l^ 

Pn  =  -  ,  (23) 

M 

where  is  the  nth  component  of  2£m*  Estimate  the  element 
gains  as 


§n 


N 


(24) 


Step  2 . 

Construct  {yi,  Y2»  •  *  •  /  using  (10)  and  (12).  Use 

for  {x,y,  flo/a}- 

Step  3 . 

Calculate 


yii!5  +  ^^  +  •••  +  ^  ^ 

M 


(25) 


Step  4 . 

Define  a  function  <P(^',£y)  with  (14)  and  (15). 

Step  5 . 

Identify  ^  as  the  which  maximizes  (P(^',Ry)  . 

The  implementation  of  Step  5  will^  be  discussed  later,  after  a 
discussion  of  the  properties  of  <P(^',^). 
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3.2  Algorithm  2:  Source  Directions  Not  Known 


In  this  method,  the  estimates  of  element  gains  are  used  to 
equalize  the  gains  in  the  array  snapshots.  These  modified  snapshots 
are  then  used  to  estimate  the  calibration  source  directions.  The 
procedure  is  as  follows: 

Step  1 . 

Step  1  of  Algorithm  1. 

Step  2A. 

Construct  a  set  of  modified  array  snapshots  {xj*  ,252  ,  •  •  •  , 2£m  }  as 
2£m'  =  S'MS)  2&n  .  m=l,2,---,M  .  (26) 

Step  2B. 

Estimate  the  calibration  source  directions  from 
{Xi  ,252  /  •  •  •  ,2&< }  •  Denote  them  by  {9^,  ,  9m}  • 

Step  2C. 

Construct  {^1,^2/ ’  ’  *  /  using  (10)  and  (12).  Use  {25n/y:m, 

for  {25,i!:/^c/a}  • 


Steps  3  to  5 

Steps  3  to  5  of  Algorithm  1. 

The  estimation  of  the  source  directions  in  Step  2B  is  highly 
dependent  on  the  array  and  the  interference  conditions.  Some 
examples  will  be  given  in  Section  IV. 

3.3  Maximization  of  (P(^',fiy) 

A  preliminary  study  of  the  properties  of  (P(<t' ,  as  a 
function  of  ,  has  been  made.  Some  notable  results  are  listed 
below: 

1.  The  element  gain  estimates  are  usually  very  accurate  when  the 
SNR  at  the  array  elements  is  30  dB  or  larger  and  the  number  of 
array  snapshots  is  reasoncQsly  large. 

2.  The  function  <P(i',|y)  is  usually  not  unimodal.  However,  if  the 
SNR  is  30  dB  or  larger,  the  global  maximum  is  usually  the  one 
closest  to 

3 .  The  finite  precision  in  computer  arithmetic  introduces  a  small 
random  error  to  the  value  of  (P(^',fiy)  calculated.  This  error 
is  very  difficult  to  handle  in  the  maximization  process  unless 
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always  satisfies  (2) .  The  error  also  generates  many  false 
maxima  in  the  neighbourhood  of  each  true  maximum.  The  density 
"and  peak  values  of  these  false  maxima  generally  increase  as 
the  Euclidean  distance  from  the  true  maximum  decreases  or  as 
the  value  of  the  true  maximum  increases. 

Because  of  the  cd^ove  observations,  Step  5  in  the  estimation 
procedure  is  implemented  in  the  following  manner: 

Step  5A. 

Initialize  the  components  of  as  samples  of  a  random  variate 
uniformly  distributed  in  the  range  [-90,90]  degrees.  Use  hill¬ 
climbing  with  random  searches  at  each  local  meocimum  in  the 
maximization.  Modify  one  component  of  at  a  time.  Impose  (2) 
on  the  modified  before  calculating  (P(^',Ry). 

Step  5B. 

Repeat  the  Step  5A  K  times,  using  a  large  K.  Denote  the 
maximum  positions  and  values  by  k=l,  2 ,  •  •  • ,  K}  , 
respectively.  Identify  the  largest  (P^  as  the  value  at  the 
global  maximum  and  denote  it  by  <P.  Identify  the  corresponding 
^  as  the  estimate 


4  PERFORMANCE  EVALUATION 


Four  examples  have  been  designed  to  evaluate  the  performance 
of  the  array  calibration  method.  The  scenarios  are  the  four 
combinations  of 


(a)  Calibration  source  directions  known 
or 

Calibration  source  directions  not  known,  and 

(b)  Number  of  interfering  sources  larger  than  number  of  array 
elements 

or 

A  strong  specular  multipath  interference  present. 

The  true  element  gains  in  the  simulations  are  unity  and  the 
components  of  ^  are  samples  of  a  Gaussian  variate  with  zero  mean 
and  eight  degrees  standard  deviation.  Five  hvindred  (K=500}  sets  of 
calculated  in  Step  SB  in  each  example. 

It  is  not  necessary  to  evaluate  the  calibration  method  with 
true  element  gains  other  than  unity.  As  long  as  the*?®  gains  are 
close  to  unity,  the  percentage  errors  in  the  imates  are 
insensitive  to  the  true  values.  This  property  can  be  verified  by 
noting  that  the  replacement  of  2=  (1, 1,  •  •  • ,  1) by  g=£,  leads  to  the 
replacement  of  2&n  by  diag{ig}3c„  and  §„  in  (24)  by  *:i8„§n.  The  scaling 
factor  K.  is  close  to  unity  if  the  components  of  ^  are  also  close  to 
unity . 

In  the  presentation  of  results,  the  maximum  relative  to  the 
largest  maximum  S*  is 

Pk  =  -  ;  (27) 

<S> 


the  root  mean  square  (RMS)  error  in  ^  is 


ev  = 


1^  -  i.\ 
Vn 


(28) 


the  distance  between  and  ^  is  ;  and  direction  in  beamwidth 

is  N  d  sin0 ,  where  d  is  the  element  spacing  in  wavelengths . 


Example  1 . 

In  this  example  the  calibration  source  directions  are  known 
and  the  number  of  interfering  sources  is  larger  than  the  number  of 
array  elements.  The  array  is  shown  in  Fig.  1.  It  has  eight  elements 
(N=8)  spaced  at  a  distance  of  half  a  wavelength  (d=0.5  wavelength)  . 
The  signal  scenario  is  shown  in  Fig.  2.  A  calibration  source  giving 
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a  30  dB  SNR  at  the  array  elements  is  used.  It  moves  in  direction 
from  -6.28  to  6.28  degrees  and  is  represented  by  a  horizontal  line 
with  arrow  heads  at  the  ends.  Ten  10-dB  SNR  interfering  signals, 
represented  by  vertical  arrows,  are  present.  Four  are  located  at  (- 
41,-19,0,19)  degrees  and  six  are  in  the  range  [33,50]  degrees.  They 
are  uncorrelated  with  each  other  and  with  the  calibration  signal. 
Four  hundred  (M=400)  array  snapshots  are  generated.  Algorithm  1  for 
known  calibration  source  directions  is  used. 


S 


Fig.  1  Array  Configuration  in  Examples  1  and  2. 


The  true  and  estimated  values  of  element  gains  and  phases  are 
shown  in  Table  1.  The  errors  in  gain  estimates  are  1.5%  or  less  and 
the  errors  in  phase  estimates  are  0.2  degree  or  less.  Out  of  the 
500  maximisation  runs  in  element  phase  estimation,  481  terminate  in 
the  region  of  the  global  maximum.  Each  produces  a  maximum  value 
between  1035.03  and  1037.01  a^,  where  is  the  white  noise  power. 
The  positions  of  the  ten  largest  maocima  are  separated  from  each 
other  and  from  the  true  ±  by  less  than  0.4  degree.  The  other  19 
runs  terminate  at  positions  with  (P^  between  9.16  and  12.01  a^.  These 
positions  are  separated  from  each  other  and  from  the  true  ^  by  at 
least  10  degrees.  The  distribution  of  relative  power  is  shown  in 
Fig.  3,  where  the  horizontal  line  between  p=0.012  and  0.998 
reflects  the  absence  of  any  (P^  between  12.01  and  1035.03  .  The 

dependence  of  Ck  on  in  the  region  of  the  global  maximum  is  shown 
in  Fig.  4.  Here,  the  RMS  error  is  less  than  0.25  degree  for  p^ 
larger  than  0.9990. 
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Table  1:  True  and  estimated  values  of  element  gains 
and  phases  in  Example  1. 
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RMS  ERROR  c  (DEG)  CUMULATIVE  FRACTION 


RELATIVE  POWER  p 


3  Distribution  of  output  powers  in  Example 
Calibration  source  directions  are  know. 


Fig.  4  Dependence  of  Ck  on  Pk  in  the  global 
maximum  region  in  Example  1.  Calibration  source 
directions  are  known. 


Example  2 

In  this  example  the  calibration  source  directions  are  knovm 
and  a  strong  specular  multipath  interference  is  present.  The  array 
is  a  radar  antenna  used  to  track  low-flying  targets  over  water.  It 
operates  at  10  GHz,  has  eight  isotropic  elements  spaced  at  a 
vertical  distamce  of  12.5  cm  (da4.17  wavelengths),  and  is  centred 
at  a  height  of  15  m.  The  arrangements  for  measuring  calibration 
data  are  shown  in  Fig.  5.  An  active  source  is  mounted  on  a  tower 
500  m  in  front  of  the  array.  It  moves  up  in  equal  steps  from  10  to 
25  m.  A  flat  earth  model  for  the  specular  multipath  is  used.  The 
reflection  coefficient  at  the  surface  is  -1,  corresponding  to  0  dB 
power  attenuation  and  180  degrees  phase  change.  The’ noise  is  white, 
and  the  SNR  at  the  array  elements  is  30  dB.  Four  hundred  array 
snapshots  are  generated.  Algorithm  1  is  used. 


SOURCE 


Fig.  5  Arrangements  for  calibration  data  measurement  in 
Example  2  and  4 . 


The  true  and  estimated  element  gains  and  phases  are  shown  in 
Table  2,  where  the  errors  in  gain  and  phase  estimates  do  not  exceed 
0.4%  and  0.3  degree,  respectively.  Out  of  the  500  maximization  runs 
in  phase  estimation,  238  terminate  in  the  region  of  global  maximum 
and  have  (P^  between  999.51  and  1000.61  a^.  The  ten  positions  with 
the  largest  maxima  are  separated  from  each  other  and  from  ^  by  less 
than  0.7  degree.  The  distribution  of  relative  power  is  shown  in 
Fig.  6  and  the  dependence  of  RMS  error  on  p^  in  the  region  of 
global  maximum  is  shown  in  Fig.  7. 


ft 
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Fig.  7  Dependence  of  on  py^  in  the  global  meocimum 
region  in  Example  3 .  Calibration  source  directions  are 
known. 

In  this  example,  the  direct  and  specular  multipath  signals  are 
equally  strong  at  the  array  elements.  The  peak  values  of  some  local 
maxima  are  larger  than  420a^.  As  the  calibration  source  moves  up 
from  10  to  25  m,  the  angular  separation  between  the  direct  and 
reflected  signals  increases  from  2.29  to  5.73  degrees  (1.52  to  3.81 
beamwidths) ,  and  their  phase  difference  at  the  array  centre 
increases  from  7020  to  17820  degrees. 

Example  3 

Here,  the  calibration  source  directions  are  not  known  and  the 
number  of  interfering  sources  is  larger  than  the  number  of  array 
elements.  The  array,  signal  scenario,  and  calibration  data  are 
identical  to  those  in  Example  1.  Algorithm  2  is  used. 

The  calculation  of  9,  in  Step  2B  is  as  follows.  Initially,  jSn! 
is  used  with  the  root  variant  of  the  modified  covariance  method 
[12]  and  a  third  order  linear  prediction  filter  to  produce  three 
candidate  source  direction  estimates.  Next,  the  power  associated 
with  each  estimate  is  calculated.  Finally,  the  estimate  with  the 
largest  power  is  identified  as  9„. 

In  the  results,  the  element  gain  estimates  are,  as  expected, 
identical  to  those  of  Example  1.  The  distribution  of  relative  power 
p,c  is  essentially  the  same  as  Fig.  3.  The  dependence  of  on  in 
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the  region  of  global  maximum  is  shovm  in  Fig.  8.  Here,  RMS  error 
is  less  than  0.75  degree  for  larger  than  0.9990.  The  ten  largest 
maxima  are  separated  from  each  other  by  less  than  0 . 3  degree .  Their 
distances  from  ^  are  slightly  less  than  2.0  degrees. 


Fig.  8  Dependence  of  e,^  on  p^  in  the  global  maximum 
region  in  Example  2.  Calibration  source  directions  are 
estimated. 


The  estimate  ^  is  given  in  Column  3  of  Table  3 .  Its  deviation 
from  the  true  ±  in  Column  2  is  the  ^  in  Column  4.  This  deviation  is 
decomposed  into  a  linear  component  ^  in  Column  5  and  a  random 
component  ^  in  Column  6,  using  the  least  squares  linear  curve 
fitting  method.  The  linear  component  decreases  in  steps  of  0.291 
degree  down  the  table .  The  rcuidora  component  is  0.2  degree  or  less 
in  magnitude. 

The  linear  component  6i,  produces  a  bias  in  signal  direction 
estimates  when  ^  is  used  in  array  calibration.  This  bias  comes  from 
the  errors  in  estimating  the  calibration  source  directions. 
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PHASE  (DEG) 

n 

6 

6r 

1 

-6.8 

-5.6 

1.2 

1.0 

2 

-1.9 

-1.3 

0.6 

0.7 

3 

7.5 

7.8 

0.4 

0.4 

-0.1 

4 

2.7 

2.9 

0.1 

0.1 

0.0 

5 

-7.1 

-7.2 

-0.1 

-0.1 

0.0 

6 

7.2 

6.5 

-0.7 

-0.2 

7 

3.3 

2.6 

-0.7 

-0.7 

8 

-4.9 

-5.8 

-0.9 

-1.0 

Table  3 :  True  element  phases  estimated  phases  errors 
in  estimates  5,  linear  component  of  errors  and  random 

component  of  errors  ^  in  Example  3 . 


Example  4 

Here,  the  calibration  source  directions  are  not  known  and  a 
strong  specular  multipath  is  present.  The  array,  signal  scenario, 
and  calibration  data  are  identical  to  those  in  Example  2 .  Algorithm 
2  is  used. 

The  calculation  of  in  Step  2B  is  as  follows.  Firstly,  the 
2£„!  is  used  to  calculate  three  candidate  direction  estimates.  Next, 
the  power  associated  with  each  estimate  is  calculated  and  the 
estimate  with  the  smallest  power  deleted.  Finally,  the  remaining 
estimate  with  a  larger  elevation  is  identified  as 

The  element  phase  estimates  are  given  in  Column  3  of  Table  4 . 
The  errors  have  a  linear  component  which  changes  by  0.535  degree 
between  adjacent  elements.  When  this  component  is  removed,  the 
error  magnitudes  are  0.3  degree  or  less.  The  distribution  of 
relative  power  is  essentially  the  same  as  Fig.  7.  The  dependence 
of  on  Pk  in  the  region  of  global  maximum  is  shown  in  Fig.  9, 
where  the  RMS  error  is  less  than  1.5  degree  for  p^  larger  than 
0.9990.  The  ten  largest  maxima  are  separated  from  each  other  by 
less  than  0 . 9  degree  and  their  distances  from  ^  are  less  than  4 . 0 
degree . 
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PHASE  (DEG) 

n 

<t> 

6 

1 

-6.8 

-5.0 

1.8 

1.9 

0.0 

2 

-1.9 

-0.9 

1.1 

1.3 

-0.3 

3 

7.5 

8.4 

0.9 

0.8 

0.1 

4 

2.7 

3.1 

0.3 

0.3 

0.0 

5 

-7.1 

-7.1 

0.0 

-0.3 

0.3 

6 

7.2 

6.6 

-0.6 

-0.8 

0.2 

.7 

3.3 

1.9 

-1.4 

-1.3 

-0.1 

8 

-4.9 

-7.0 

-2.1 

-1.9 

-0.2 

TcQsle  4 :  True  element  phases  estimated  phases  ^ ,  errors 
in  estimates  i.,  linear  component  of  errors  and  random 

component  of  errors  ^  in  Example  4 . 


RELATIVE  PONER  p 


Fig.  9  Dependence  of  on  in  the  global  maximum 
region  in  Example  4 .  Calibration  source  directions  are 
estimated. 
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5  CONCLUSIONS 


A  method  which  uses  a  movable  calibration  source  to  estimate 
the  element  gains  and  phases  of  a  sensor  array  has  been  described. 

This  method  is  highly  tolerant  to  interference,  including  specular 
multipath.  The  number  of  interfering  sources  also  need  not  be  known  ^ 

or  be  smaller  thcui  the  number  of  array  elements .  Computer 
simulations  with  reasonably  strong  calibration  signals  produce 
highly  accurate  gain  and  phase  estimates  when  the  source  directions 
are  known.  When  the  directions  are  not  known,  the  set  of  element 
phase  estimates  contains  an  additional  component  which  corresponds 
to  a  bias  in  direction  estimates. 

This  method  has  many  applications.  One  example  is  the 
calibration  of  a  high-frequency  array  used  in  direction-finding. 

The  method  accounts  for  the  presence  of  other  signal  sources  as 
well  as  the  indirect  signals  from  the  ionosphere  and  other 
reflectors.  Another  example  is  the  calibration  of  a  sonar  array  in 
the  ocean.  The  method  accounts  for  the  surface  and  bottom 
reflections  as  well  as  the  interfering  signals  from  ships  and 
marine  life.  A  third  example  is  the  calibration  of  an  array  antenna 
used  to  track  low-flying  targets  over  water.  The  method  accounts 
for  the  specular  reflection  at  the  water  surface.  ^ 

The  calibration  method  has  the  limitation  that  mutual  coupling 
among  the  array  elements  is  absent.  It  appears  that  there  is  no 
simple  modification  which  can  take  this  coupling  into  j 

consideration . 
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of  element  phase  estimates  contains  an  additional  component  which  produces 
a  bias  In  direction  estimates. 
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